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Abstract 

We have adapted R-matrix theory to calculate phonon scattering across systems of molecular 
to mesoscopic scale. The key novelty of this work is that the only required information about 
the scattering region are its normal modes, which are evaluated only once for a system. Thus, 
R-matrix theory is a computationally efficient and simple approach to calculate phonon scattering 
in larger systems. To validate and to demonstrate the applicability of the theory, we apply it to 
two systems: a one dimensional chain of atoms and a graphene nanoribbon. In both cases, we 
discuss the effect of mass impurities on thermal transport. 
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NOMENCLATURE 



a- lattice constant 
A- constant 

A- matrix defined for convenience 
b- basis index 

B- matrix defined for convenience 
c- phonon velocity 
h- Planck constant 

h- reduced Planck constant (h = h/2n) 

IR- interior region 

X- identity matrix 

%- unitary imaginary number 

j- lattice index 

/C- mass normalized force constant matrix 
ks- Boltzmann constant 
/- lead index(l,2) 
L- Lead 

Cb- Bloch operator 
m- mass 

mira- Impurity mass 

Ai- dynamical matrix 

n- integer use to number the eigen modes 

N uc - number of atoms in a unit cell 

N^- number of phonon branches at frequency 00 

p- phonon branch 

q- wave vector 

Q- thermal current 

1Z- R-matrix 

S- scattering matrix 

t- time 
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T- temperature 

u- displacement of an atom 

u- displacement vector 

v- eigen vector in the interior region 

V- many coordinate inter-atomic potential 

x- dimensionless variable 

Greek Letters 

a,/3- Cartesian degree of freedom 
r~ transmission probability 

e- polarization vector of atoms in a unit cell in the lead 

r\- Bose-Einstein distribution function 

cu- phonon frequency 

Co- dimensionless phonon frequency 

ojd- cut off frequency 

a- thermal conductance 

a- dimensionless thermal conductance 

A- eigen value in the interior region 

t- dimensionless temperature 

0- force constant 

X~ column vector defined for convenience 
£- unit vector in 3N UC dimensional space 

Superscripts 

* - complex conjugate 
f - Conjugate transpose 
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I. INTRODUCTION 



The rapid advancement of nanotechnology drives the fabrication of structures much 
smaller than the mean free path of electrons as well as phonons. Electronic transport on 
the nanoscale has been studied for over three decades and fascinating quantum effects have 
been observed^-. Phonon transport on this scale is of significant interest because of the 
increased power dissipation in nanoelectronics, which undermines the correct functionality 
of devices and limit their life time^ - -. The need for better thermal conductivity has driven 
the search for a better understanding of thermal transport on the nanometer level. 

For example, the existence of a quantum of thermal conductance in mesoscopic systems 
has been proposed 9 and demonstrated experimentally^ 1 ^. In nanowires (width < 100 nm), 
phonon transport is one dimensional at low temperatures due to the lateral confinement.— 
In the ballistic limit, where the wire length is less than or equal to phonon mean free path, 
thermal conductance reaches to its maximum ir 2 k 2 B T /?>h, which is called the quantum ther- 
mal conductance. In this limit, thermal transport is largely determined by the boundaries 
(where the boundaries are defined by the geometry of the device), interfaces and impurities. 
Thus, it is essential to calculate the thermal conductivity including all the details of the 
atomic constitutes to include the above effects. 

Analogous nanoelectronic transport has been studied using the Landauer and Buttiker 
formalism, in which the transport is determined by the transmission probabilities^"- 1 ^^. 
The phonon transport in these systems can also be treated by the very similar approach 9 - 1 ^. 
In this context, the theories that was using initially in calculating electronic transmission 
probabilities should be applicable to phonons. For instance, the Green's function approach 
has been successfully developed to calculate phonon transmission probabilities^ - —. 

R-matrix theory is a theoretical approach commonly used in nuclear and atomic physics 
to solve scattering problems^ - —. The first development of the R-matrix theory draws back 
to the work of Wigner and Eisenbud in nuclear scattering in 1947. Recently, this approach 
has been successfully developed to calculate electronic scattering in mesoscopic quantum 
devices^ - —. In this work, we present the first mathematical development of the R-matrix 
theory to calculate phonon transmission probabilities across a device of molecular to meso- 
scopic scale. 
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The key novelty of R-matrix theory for phonon is that the only required information about 
the scattering region is its normal modes, which is needed to evaluate only once for a system. 
Normal modes are the collective oscillations of the lattice, which are orthonormal and form a 
complete basis in a finite region. These are calculated by diagonalising the dynamical matrix 
in the harmonic approximation of the interaction potential. The harmonic approximation 
is well suited to study thermal transport in mesoscopic scale systems and low temperature 
transport, which is not accessible by the classical molecular dynamic simulations. R-matrix 
theory is an efficient and simple way in calculating the phonon transmission through a 
nanoscale system. 

This paper is structured as follows: Section [TT] describes the mathematical formalism of 
the R-matrix theory. An introduction to the model system and the scattering picture of 
the lead is in section III. II In section III.2| we construct the phonon version of the R-matrix. 
Then we solve for the scattering matrix (S) in section III. 31 To validate the theory and to 
demonstrate the applicability, we apply the theory for two type of systems (Sec. IIIII) . First, 
we apply the theory to a one dimensional chain of atoms. Then we calculate the phonon 
transmission through a graphene nanoribbon. In both cases we discuss the effect of mass 
impurities to the thermal transport and a conclusion is given in section IIVI 

II. THE R-MATRIX FORMALISM 

II. 1. The system and the scattering picture 

The essence of R-matrix theory is to divide the system into an interior regioniJR) and 
two or more asymptotic regions. The IR is the main scattering center of the system. The 
asymptotic regions are the semi-infinite leads connected to the IR (FigJT]). In the asymptotic 
regions, we assume the leads are narrow enough that the phonon transport in leads is strictly 
along the length direction of the leads. In the lateral direction there are confined phonon 
modes which is also called the phonon subbands. Thus, we define a extended unit cell, 
depicted by the boxed regions in the figiH that has the periodicity in the length direction. 
Although R-matrix theory can handle any number of leads, for simplicity we consider the 
case of only two leads, while it is straight forward to generalize to a system with many leads. 
The IR and the leads are seperated by a set of boundaries in the periodic region, which 
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FIG. 1. Schematic of the model system having two semi- infinite leads connected to a device. The 
leads are assumed to be periodic and two boundaries, B\ and B2, are defined in the periodic region. 
The region between boundaries is the interior region and the region outside is called the asymptotic 
region. The unit cells in the leads are depicted by the boxed regions. The lead lattice points are 
indexed symmetrically by ji, where I = 1(2) refers to the lead 1(2). 

we denote by B\. The unit cells in the leads are indexed by ji, where Z = 1(2) is the lead 
1(2). The lattice point ji = is at the boundary Bi in to the IR and ji increases positively 
outward the IR. The heat current from lead 1 to 2 (^1,2) is given by the Landauer Formula 

Qi,2= —hu(r] Ll {uj)-T] L2 {uj))T lfl (uj), (1) 

/, , ^ 

" < J -'m,in 

where, 77l ; (w) = {exp (hu/kBTi) — l} -1 is the Bose-Einstein distribution function of lead I, 
ks is the Boltzmann constant, h is the reduced Planck constant and T] is the temperature 
of lead I. The transmission probability of a phonon of angular frequency u from lead 1 to 2 
is given by Fi^uj). 

Since the leads are periodic there are well-defined phonons in each lead. The details of 
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the phonons can be obtained solving the dynamical equation, 



a,/3 p 



(2) 



b'/3 

over a unit cell, where the dynamical matrix M. is, 




(3) 



In the harmonic limit, force constants can be calculated by, 



(4) 



which is evaluated at the equilibrium positions of atoms. The symbol V refers to the many- 
coordinate inter-atomic potential. The coordinate of the atom b in the unit cell ji is given 
by Rj,,6 = Rj; + R&, where R^ is the position of the lattice point ji and H b is the position 
of the atom b within a unit cell. The Cartesian degree of freedom is represented by a and 
(3, m b is the mass of the atom b, q is the wave vector, m" is the displacement of the atom j 
at the direction a, and e£ is the polarization amplitude of atom b in direction a in phonon 
branch p. There are 3iV uc phonon branches in each lead, where N uc is the number of atoms 
in the unit cell. 

We can construct a scattering solution in the leads with these phonon modes. The 
scattering solution is comprised of three components: the incoming wave, reflected wave 
and transmitted wave. We consider an incoming phonon of frequency u, whose energy 
is hu, in lead Iq and branch p . It can reflect in to any phonon branch (p) available at 
u in lead l and transmit in to any phonon branch (p) in lead /, which is different from 
l . We assume all the scattering processes are elastic and the anharmonic interactions are 
neglected. As we defined our leads symmerically, the incoming waves are going at the 
direction of decreasing the lead lattice index ji and the scattered waves (both reflected and 
transmitted) are going at the increasing direction of the lattice index. Thus, the incident and 
scattered components of the scattering solution can be written in terms of plane waves as, 
q p (cu)) and X p (q p (^)) respectively, where X p (q p (w)) is a vector whose entries are the 
values of e£exp(— iut + iaf (uj) Jlj ub ) / ^/m^ of all the atoms in a unit cell. Now the scattering 
solution (uj o ) of the lead / for above described scattering event can be written as, 





p 
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where, the first part inside the summation represents an incoming wave in branch po in lead 
lo and the second part is the scattered wave. The scattering amplitude in branch p lead I is 
given by S 1 ^ (oj), which is a complex number. The scattering amplitudes form a 2iV w x 2N UJ 
scattering matrix S(u), where the N w is the number of phonon branches at frequency (u). 
The matrix elements of S, when I = l are the reflection coefficients and / 7^ Iq are the 
transmission coefficients. 

From the scattering amplitudes, we can calculate the transmission probabilities. The 
transmission probability is the ratio of the energy fluxes carried by each vibrational wave. 
The transmission probability from branch p in lead Zq to branch p in lead Iq) can be 



expressed as, T 



i P 

I0P0 



,00 



S\fl ( w )SifpS u ^ cP ( w } I cP ° where c p (u) is the phonon velocity 



of the branch p at frequency uj— . Thus the transmission probability from lead Iq to Iq) 



is. 



(6) 



II. 2. R-matrix 



The dynamics of the entire system can be described by an infinite set of coupled equations, 
which are usually written as matrix form, 



(w 2 X + /C) u = 0, 



(?) 



where X is the identity matrix, u is the displacement vector and /C is the mass normalized 
force constant matrix. The matrix elements of /C can be calculated as, 

1 



id' 



IJljlTlji 



(8) 



where j refer to a common index of atoms used to label all the atoms in the system. The 
matrix /C can be expressed more descriptively using the block matrices: Kir is defined only 
in the finite IR, IC^ is defined in the lead I that is infinite, and K IR Ll is the coupling between 
the interior region and the leads. 



/ 



/Cu 



KlR,L! fclR fclR,L 2 
K-L 2 ,IR K-L 2 J 



(9) 
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where and Ui R are the displacement vectors of atoms in lead / and the IR. 

It is natural to describe the IR by normal modes of vibrations. However, in our scattering 
problem the IR is connected to semi-infinite leads and thus we can not solve for the normal 
modes. The solution is to solve the equations of motion for a different, but related physical 
system. When solving the continuum Schrodinger equation in a finite region, a "Bloch oper- 
ator" is included to solve the problem of non-Hermiticity of the kinetic energy operator— 
The resulting Hamiltonian is called the "Bloch Hamiltonian" , which is explicitly Hermitian 
inside a finite region. In the phonon problem we define a "Bloch operator (£#)" , which 
includes the couplings between the IR and the leads as follows, 



K, LuIR 

K>L 2 ,IR 



(10) 



The term C^u gives the coupling forces between the IR and the leads. It can involve 
interactions of any finite range, e.g. nearest neighbor, next nearest neighbor, etc. We 
subtract the term from the both side of the eqJT] yielding, 



(cu 2 I + fC — Cb) u = -£gu. 



(11) 



The matrix (JC — Cb) on the left is called the "Bloch dynamical matrix", in which the 
coupling between the two regions are removed. The eq{TT] can be written in block matrices 
as follows, 



ULi 



K T « o 



\ 



V 



K IR 
K 



Li J 



K. 






K>Li,IR 





IRM 





fclR,L 2 





K>L 2 ,IR 






w„ \ 



• (12) 



We can extract the center row of the eqfl2l which describe the IR, 

(u 2 u IR + K,irU JR ) ^ C B l u Ll , 



(13) 



where Cg 1 = ICj Ri i r In the "Bloch dynamical matrix" IR can be solved independently. We 
find eigenvalues (A n ) and eigenvectors (v( n )) of JCi R according to, 

(X 2 n + K IR )v^=0. (14) 
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These eigenvectors are orthonormal and form a complete set as the mass normalized force 
constant matrix is always Hermitian. A general solution in the IR can be expanded using 
the above eigen vectors as follows, 

u IR = J2^ {n) - (15) 
By plugging this in to the left hand side of the equation (!T3|) yields, 

J2(^-^ 2 )A n ^ = J2^ l u Lr (16) 

n I 

By taking the inner product with y( n )t (t refers to the conjugate transpose) we find, 

I X n~ U 

Finally, the scattering solution in the IR can be expressed as, 

u ^ = E v(n) X 2_ U 2 ■ ( 18 ) 



n,l 



it 



According to the eqJH] and \TU\ Bloch operator of lead I, Cq 1 , has infinite number of 
columns and 3iV rows, where N is the number of atoms in the IR. However there are non- 
zero couplings only between a subset of atoms in the IR with a subset of atoms in the leads, 
in the visinity of the boundary B[. Since the boundary is defined in the periodic region, we 
assume that the minimum size of the coupling regions to be unit cells at either side of the 
boundary. This is a reasonable (but not essential) assumption for most practical situations. 
Therefore we can construct a matrix of dimensions 3N UC x 3N UC that contains the non-zero 
couplings in Lb . This effective part of the Bloch operator is denoted by Cb ■ Now we replace 
by uj (ji = 1), of which the entries are the displacements of the atoms in the unit cell 
at lattice point ji = 1. Similarly (~v^) l contain the normal modes amplitudes of the atoms 
in the unit cell at ji = 0. Now we can replace v n ^ (Cb^Li) by (v^ n ^) ; (/lj3 l u l lopo (ji = 1) \ in 
eqUHJ 

In the electron transport problem, the R-matrix relates the value of the wave function to 
its normal derivative at the boundary of the scattering region and the normal derivative at 
the boundary serves as the Bloch operator.— The phonon version of the R-matrix (1Z) can 
be defined as follows, 

* v 'Ui» = o) = E rX/r'^V V(j, = i), (19) 

s",s',l' 
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n 11 

where £ s, s are a set of unit vectors that are ortho normal and complete in the 3iV nc di- 
mensional space. These vectors can be chosen arbitrarily and a possible choice is given in 
appendix lAl 

The eqJUj] is obtained by matching the lead solution and the interior region solution at 
the unit cell at ji = 0. In general, the matching is done on the surface atoms that share 
the both regions. The choice of the unit cell is adapted for the convenience of explanation. 
This particular choice is valid for some practical cases such as zigzag graphene nanoribbons. 
Moreover, the incorporation of unit vectors £ s is not essential. It allows to calculate the 
matrix elements of the 1Z before hand so that it can be used in the scattering matrix 
calculation. Alternatively, one can keep the summation in eqJT9]only over on "n" and 
without including the unit vectors £ s and proceed. 



II. 3. Scattering Matrix 

The scattering matrix (S) can be obtained by plugging the explicit form of the scattering 
solution in the lead (u| opo ) from eqj5]in to the R- matrix equation ( [191) . We can express the 
scattering matrix as, 

S{u) = - H(q) - B(ci)]- 1 . lA(-q) - B(-q)] , (21) 
where the matrix elements of A and B are given by, 

Equation (121 p is our key expression for the scattering matrix. A derivation of this is given 
in appendix [Bj The matrices A and B are independent of the number of atoms in the IR. 
It is the R-matrix that represent the IR in the scattering calculation. In order to construct 
the matrices A, B and 71, we only need to calculate for the normal modes of the IR in the 
Bloch dynamical matrix and the lead phonons. The appealing feature is that these details 
are needed to be calculated only once for a system. 
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III. APPLICATIONS 



III.l. One dimensional chain of atoms 

First, we calculated the transmission probabilities for a one dimensional germanium chan- 
nel connected to silicon leads using R-matrix theory and could reproduce the results in 
reference 1 ^ which was calculated using the nonequilibrium Green's function method. Here, 
we use R-matrix theory to calculate the phonon transmission of a one-dimensional (ID) chain 
of atoms that has one atom in the unit cell and investigate the effect of mass impurities to 
the thermal transport. In particular, we study three cases: a single mass impurity, diatomic 
chain, and a random mass impurity distribution connected to two semi infinite uniform 
leads. 

The atomic masses are rrij = 1.0 for all j and the couplig constants are given by, 



The IR of the system consists of 25 atoms. The Bloch operator (C l B ) for this ID system 
is a single number, which is the force constant between the atom in the boundary and the 
adjacent atom in the lead. We consider only the longitudinal vibrational modes. 

The dotted line in figfS^, shows the transmission probability for uniform chain. It gives 
unit transmission for all the frequencies up to the cut off up . Mass impurity is made by 
changing the mass of an atom in the middle of the interior region to mfl = 2.0. The 
dashed line in figj2^, shows the transmission probability calculated in the presence of a single 
mass impurity. In the semi-classical approach impurity point defect scattering is modeled 
by a Rayleigh type expressionism that is proportional to w 4 . We observe a similar trend 
in the transmission probability for scattering from a single mass impurity in ID. Then we 
incorporate another mass of rnf^ = 0.6 with the forming a periodic diatomic structure 
in the middle of the interior region. The number of unit cells included is seven. The 
corresponding transmission is plotted by the solid line in figfj^i. We can see the opening 
of a forbidden region having zero transmission. Next the above two masses are randomly 
distributed replacing 14 atoms in the interior region. In this case the transmission is depicted 



2 for 



1 — 3 




< — 1 for j = i ± 1 



(23) 







otherwise 
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FIG. 2. (a) Transmission probability (T) as a function of dimensionless frequency (u/wd). (b) 
Dimensionless thermal conductance (a) as a function of dimensionless temperature (r). The dotted 
line is for the uniform chain. Dashed line is for the inclusion of single impurity of rri^l = 2.0 and 
periodic diatomic structure of 7 unit cells with masses and ttlj^L = 0.6 is given by the solid 
line. The thick dotted dashed line is for the random distribution of 14 atoms with masses 
and mf m . 



by the dotted dashed line. 

Now we calculate the thermal conductance as a function of temperature. In the linear 
response regime, where AT < T, AT = \T X - T 2 | and T = {T 1 +T 2 ) /2, the thermal 
conductance (a) of the system can be expressed as follows, 



klT [ Xma * , F(x)x 2 e x 
IT I ' 

** <A* 771.7. 



(e x - 1) 



2 ' 



(24) 



where, x = -fi^. The dimensionless thermal conductance (er = cr/(kBU)D)) is denned as, 



1 f 1 , T{Cj)Cj 2 e*l T 
a — — I doj- 



(25) 



2vrJ "V^/t - l) 2 ' 

where r = (hsT/hup) , ojd is the cut off frequency, which is equal to 2.0 in this problem , 
Cj = uj/ujd- 
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The dimensionless thermal conductance for previously described cases is plotted in figj2b. 
The random distribution of masses shows the lowest thermal conductance. The random 
distributions could creates localized states suppressing the phonon transport whereas the 
periodic structures are promising for travelling waves. 

III. 2. Graphene nanoribbons 

Graphene nanoribbons (GNRs) are 10 ~ 100 nm wide strips of graphene with high aspect 
ratios. Graphene is a promising material for nanoscale applications due to its exceptional 
electronic^ and thermal properties^ - —. Thermal conductivity of graphene is found to be 
greater than 3000 W/mK in recent experiments^ - —. Moreover it has been recently revealed, 
both theoretically^ and experimentally^, that the flexural (out of plane) vibrational modes 
are the governing energy carriers in graphene. We apply R-matrix theory to calculate the 
transmission of flexural vibrational modes of a GNR. 

Figure [3] shows a schematic of a zigzag graphene nanoribbon (ZGNR) of six zigzag chains 
in width. A unit cell is depicted by the boxed region and the lattice constant is given by 
a. The force constants are calculated using a second generation Brenner potential^. When 
considering interactions up to the next nearest neighbors, there are non-zero couplings only 
between the atoms in adjacent unit cells. This support the assumption we make in defining 
the effective part of the Bloch operator in sec JII.2[ First we construct the K, matrix as in 
eqJH] including few unit cells in the leads. Then we extract the ICm and C l B as describe in 
the eqUanddOl 

First, we calculate the phonon dispersion of flexural modes according to the eqj2j Figure 
shows the phonon dispersion along the direction of the ribbon length, where q u is the 
wave vector projected on to the length direction. There are 12 phonon branches available 
due to the fact that the unit cell has 12 atoms. These phonon branches represent different 
laterally confined phonon subbands that are specially available in nanowires. The acoustic 
branch is also available in large graphene sheets. The lowest lying optical branch, which has 
a twisting character (mode 2 in figE]), can only exist in nanoribbons. In larger graphene 
sheets twisting modes do not present since it is energetically costly. Figure Hb shows the 
calculated transmission probability (T(u)) from the R-matrix theory. The T counts the 
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FIG. 3. Schematic of the a zigzag graphene nanoribbon (ZGNR) of six zigzag chains in width. A 
unit cell is depicted by the boxed region. The lattice constant in the leads is denoted by "a" . The 
horizontal arrow shows the direction of q n . 

number of phonon branches available at each frequency. Similar behavior of the transmission 
probability has been reported by the non-equilibrium Green's function calculations. 

We incorporate a larger mass (mass impurity of m/ m = lOOm^) replacing one carbon 
atom in the center of the interior region. We keep about 20 unit cells (240 atoms) in the 
interior region. Although this is not a realistic situation, could be analogous to a strong 
coupling to a heavy molecule. The calculated transmission probability for this case is shown 
by the dashed line in figSt 1 - We observe that the overall transmission of phonons is reduced 
due to the heavy mass. In figfS] we plot the thermal conductance (a) as a function of the 
temperature. Thermal conductance of the ZGNR with the mass impurity (dashed line) is 
considerably lower than that of uniform ZGNR (solid line). 

In R-matrix theory we are working on the real space displacement waves. This allows us to 
easily sketch the profile of the displacement waves during the scattering. This visualization 
is important to gain clear insight on to the scattering process and tuning the mode specific 
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FIG. 4. Phonon dispersion (a) of fiexural vibrational modes of the ZGNR. The transmission 
probability (T) calculated for the perfect ZGNR from R-matrix theory (b). The V count the 
number of phonon branches available at each frequency. The dashed line refers to the case where 
a larger mass in the middle of the IR. The frequency is measured in the units of cm" 1 . 



transmission of phonons. In figd we sketch the displacement wave in the three regions: lead 
1, interior region and lead 2, at frequency 50 cm -1 when the mass impurity is in the center 
of the interior region. We use the eqf5] for the leads and the eq[l8] for the interior region. 
The displacement waves in three regions matches nicely at the boundaries. There are two 
phonon subbands available at this frequency as shown in figEl The sketch in figJTl is for 
the waves coming in mode 1, which is the acoustic mode. We observe that the propagation 
of the acoustic mode is considerably attenuated with the presence of the mass impurity. In 
the acoustic mode, all the atoms in the unit cell vibrate in-phase. Thus, the mass impurity 
prevent this in-phase motion. However, we observe that the twisting mode is not affected 
by the presence of the mass impurity in the middle as implied by the fig|7]D. This could be 
attributed to the asymmetric polarization profile of the twisting mode (mode 2 in figJE]). 
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FIG. 5. Thermal conductance due to out of plane vibrations of the ZGNR as a function of tem- 
perature. Solid line is for the perfect structure and the dashed line is the result due to the larger 
mass in the center. 

IV. CONCLUSION 

We have adapted the R-matrix theory to calculate phonon scattering on the atomic level. 
The major advantage of the R-matrix approach is the reduced computational time. Moreover 
the ability to express the scattering matrix by a simple analytical expression is also an 
advantage of this approach as is the capability of easily plotting the scattering displacement 
wave. The application to devices with three or more leads is straightforward and only 
involves increasing the range of the lead index. We discuss results applying R-matrix theory 
to a ID atomic chain and a graphene nanoribbon. We observe a lower thermal conductivity 
for a random distribution of masses than the periodic arrangement in a ID atomic chain. For 
a single mass impurity, the profile of the transmission probability shows a trend similar to 
the Rayleigh type scattering. However, phonon transmission through GNRs in the presence 
of a mass impurity shows strong dependence on the mode profile of the laterally confined 
vibrational modes. 

Furthermore, there are a number of important differences between R-matrix theory for 
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FIG. 6. Two lowest lying flexural phonon branches of the ZGNR. The fig. (a) shows the phonon 
dispersion in the low frequency region. Figure (b) shows the polarization profile of atoms in an 
unit cell calculated at 50 cm -1 of the acoustic branch (mode 1) and the twisting branch (mode 2). 

electrons and for phonons. First, because the electron case is defined over a continuous 
spatial degree of freedom, there are an infinite number of eigenstates of the Bloch Hamil- 
tonian. This means that any solution involves a high energy cut-off given by the largest 
eigenvalue considered in the Bloch Hamiltonian. In the phonon case, there is no sense to 
oscillations that have a wavelength shorter than the distance between atoms in a single unit 
cell. Therefore, there are no more than 3iV modes associated with the interior region, and 
such convergence is not an issue. The electron case also suffers convergence issues related 
to the choice of basis used for solutions of the Bloch Hamiltonian. Since the scattering 
solution can have any amplitude and derivative on the boundary, it is important that the 
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FIG. 7. Displacement wave in three regions (lead 1, interior region and lead 2) is sketched together 
when a larger mass is present in the center of the interior region. The left boundary of the interior 
region (Bi) starts from the origin of the longitudinal coordinate (marked in units of A) and the 
right boundary (B2) is at 20 unit cells away from the B\ (at 40 A). The heavy mass is at the lO*' 1 
unit cell from the B\ (at 20 A). The sketch at top (a) is for waves coming in mode 1 at u = 50 
cm . The sketch in bottom (b) is for the waves coming in mode 2 at the same frequency. 

basis functions share the same feature. This motivates a variety of different approaches to 
choosing a basis, such as variational basis functions. In contrast, in the phonon case we 
merely solve for the normal modes of the system, and there is no such ambiguity. 

The electron problem has the advantage of well-defined boundaries and a clear potential 
energy. Setting up the Bloch dynamical matrix requires knowledge of the force constants 
(eq.flU)) evaluated at the equilibrium position of the atom. This may require simulating 
the system and allowing it to relax, minimizing the internal strain. Similarly, to define the 
"lead solutions" we must assume that the lead lattices are semi-infinite and uniform. In real 
systems the interior region may exert a strain that distorts the leads; the interior region 
must be large enough so that plane wave solutions in the leads are valid. 

The electronic scattering problem involves many "channels." In atomic physics these 
are usually different angular momentum states for the scattering electron, while in the 
mesoscopic device case they correspond to different subbands in a lead. Such inter-subband 
scattering is present in phonons as well. Electronic interband transitions require a periodic 
lattice potential in the leads and/or the interior region, so that it is meaningful to discuss 
electron bands. The phonon system may also have different bands if the unit cell has a multi- 
atom basis. This opens up the possibility of phonon transduction, where optical phonons 
can be turned into acoustic phonons and vice versa, through the suitable design of an interior 
region. 
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Appendix A: A choice of vectors £ s 

The set of vectors £ s can be chosen arbitrarily. One possible choice is, 
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Appendix B: Solving for scattering matrix (5) 



In this appendix, we discuss the algebra of solving for the scattering matrix (5). By plugging 
u' opo (j/ = 0) and U-\ (ji = 1) from eqj5]to eq IT9l and taking the projection of both side on to e p , 



we can obtain, 
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We define a matrix (B) for the convinent as follows, 

= E ^'V^' (^4tfV;H*=i)) (B2) 

s"s' 

We rewrite eq JB 1 1 rearranging the parts and inserting B. 
Y / e^x P ''\-q p, ''^j l =oK>> Sm -Y^ 

p'" i'p' i' P > p>» 

(B3) 



Let's focus on the R.H.S. of the eq JB3l It can be written as 

' -<We p 'V(-< 
(X(q)-B(q))5, 



'opo / 



(B4) 
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where A\y = e p "^X p '(q p '; ~R-j t =o)fii"i'- Similarly, we can show the L.H.S. of eq lB3l is 
(A(— q) — B(— q)). Thus, we can finally write, 

S = - [A(q) - S(q)]- 1 . LA(-q) - B(-q)] (B5) 
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